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ABSTRACT 

Let  .*}  be  a  set  of  n  x  n  complex  matrices  A  which  satisfy  the 
condition  I!  (I  -  zA)  1 1I  <  K/(1  -  !z|)a+1  for  some  a  >  0  and  all  |z|  <  1 

Then  it  is  shown  here  that  there  exists  a  constant  p(a,n)  such  that 

v  a 

II A  II  <  Kp(a,n)v  ,  u  =  0,1,...  .  This  forms  a  generalization  of  the  Kreiss 
resolvent  condition  (a  =  0). 


AMS  (MOS)  Subject  Classifications:  15A60,  39A1 1 

Key  Words:  a-Stable  set 

Work  Unit  Number  1  (Applied  Analysis) 


y:  f-c-rei 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-R0-C-0041 . 


SIGNIFICANCE  AND  EXPLANATION 


Let  A  be  an  n  x  n  complex  valued  matrix.  A  standard  and  useful 
result  in  matrix  theory  claims  that  all  powers  of  A  are  bounded  if  and  only 
if  the  spectral  radius-p(A)  -  is  less  or  equal  to  one  and  for  all  eigenvalues 
A  of  A  such  that  |X|  =1  the  matrix  (I  -  z A)-1  has  a  simple  pole  at 
z  =  If  we  consider  a  more  general  problem,  namely  when  the  powers  A, 
v  <■  0,1,...,  grow  at  most  as  va,  where  a  is  a  positive  integer,  then 
this  condition  holds  if  and  only  if  p(A)  <  1,  and  for  all  eigenvalues  X 
of  A  such  that  |X|  =1  the  matrix  (I  -  zA)-1  has  at  most  a  pole  of  order 
a  +  1  at  z  =  X. 


In  the  early  sixties  H.  0.  Kreiss,  while  studying  stability  of  numerical 


schemes  for  partial  differential  equations,  considered  a  generalization  of  -the 
Ox 

first  problem, described  ahove.  Namely,  given  a  set  A  of  n  x  n  complex 

' \  z  \  '  <’ 

valued  matrices,  when  all  powers  of  aV  A  are  uniformly  bounded.  These  sets 
-  called  the  stable  sets  -  were  completely  characterized  by  Kreiss  by  giving 
three  equivalent  conditions.  '■>  i  ,  \  ,  \  -V  h  ,  , 

'  -  t  .  1 

In  this  paper  we. consider  -a-stable  sets  /4(C*  i  0),  such  that  .for  any 


A  ('A  the  powers  are  uniformly  bounded  by 

'  o  .  V 


X  CAWVxa 

.  We  generalize  the 


generalis 

O'  v  h 

Kreiss  resolvent  condition  for ''^-stable  sets.  It  seems  that  ^-stable  sets 


are  related  to  the  concept  of  weakly  stable  numerical  schemes  for  partial 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MFC,  and  not  with  the  author  of  this  report. 


A  GENERALIZATION  OF  THE  KREISS  MATRIX  THEOREM 


Shmuel  Friedland 


1 .  Introduction 

In  various  instances  one  deals  with  iterative  systems  of  equations 

(1.1)  x<i+1)  =  Ax(i>,  i  =  0,1,2,...  . 

Here  e  cn,  A  c  Mn,  where  Cn  is  the  set  of  n  column  complex  vectors  and  Mn  is 

the  set  of  n  x  n  complex  matrices.  Clearly 

(1.2)  x<i>  =  aM°> 

and  thus  in  order  to  investigate  the  behaviour  of  x^ ^  for  large  i  one  needs  to  studv 
the  powers  A1',  i  «  0,1,...  .  Let  A  be  a  set  of  n  x  n  matrices.  A  is  called 
an  a-stable  set  if 

(1.3)  IIAVII  <  Kv**,  V  =  0,1,2,...  . 

Here  a  is  a  nonnegative  number  and  11*1!  is  a  norm  on  M_.  The  concept  of  stability  of 
the  numerical  schemes  for  solutions  of  partial  differential  equations  is  intimately 
connected  with  the  notion  of  stable  sets.  Consult  for  example  Kreiss  [1962) ,  Richtmyer  and 
Morton  [1967)  and  others.  It  seems  that  a-stable  sets  are  related  to  the  concept  of  weakly 
stable  numerical  schemes  for  partial  differential  equations.  See  Kreiss  [1962)  and  Forsyte 
and  Wasow  [I960].  The  stable  sets  were  characterized  completely  by  Kreiss  [1962].  In  this 
paper  we  generalize  the  Kreiss  result  to  a-stable  sets. 

Theorem  1 .  Let  a  be  a  nonneqative  number  and  A  be  a  set  of  n  x  n  complex  valued 
matrices.  Then  the  following  two  conditions  are  equivalent. 

(A)  There  exists  a  constant  K(>  1)  such  that  for  all  *  !  i  (1.3)  holds . 

( R )  There  exists  a  constant  K ( >  1 )  such  that  for  all  A  s  A 

(1.41  H(l  -  zA)*1!!  <  K(1  -  )z|)*(a+1i,  | z I  <  1  . 

The  implication  (A)  (R)  is  obvious.  The  implication  (R)  ==>  (A)  is  a  consequence 

of  Theorem  2  which  estimates  the  Maclaurin  coefficients  of  a  certain  family  of  rational 
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we  get 

(2.8) 

Note  that 


1 a  v  I  < 


[  max  |  f(z)  |]  r  ^  <  K(1 
lzl»r 


.  -a  -v 
r)  r 


min  ( 1 
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r)“V 
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,  -a  -v, 
r)  r 
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fi  ♦  S) 


This  establishes  the  first  inequality  in  (2.6).  To  obtain  the  second  inequality  in  (2.6) 

1  u 

choose  in  (2.8)  r  ■  — - — -  and  use  the  well  known  fact  that  ; 1  +  — 1  <  e.  To  see  that 

V  +  I  V ' 

(2.6)  is  sharp  for  each  v  consider  the  polynomial 


(2.9) 

Let  B  be 
operator. 
(I  -  zA)'1 


p(  z)  * 

a  Banach  space  with  a  norm 
Suppose  that  the  spectrum  of 
in  power  series 


v  a  ‘ 

«•».  Assume  that  A  s  B  *  B 
A  lies  in  the  unit  disc. 


□ 

is  a  bounded  linear 
Then  expanding 


•1  r  V  V 

(2.10)  (I  -  zA)  «  [  z  A 

v-0 

we  get 

(2.11)  Au  »  (2*i)-1  /  (I  -  ZA)-Idz  . 

I z | *r< 1 

Thus  if 

(2.12)  MI  -  zA)'1#  <  K(1  -  lz|)"a,  I  z  I  <  1 

applying  the  results  of  Lemma  1  we  obtain 

(2.13)  II  A11  !l  <  Ke(v  ♦  1  )  a  . 

It  is  an  open  problem  whether  the  estimate  (2.13)  is  sharp  in  some  infinite 
dimensional  Banach  space.  The  following  result  enables  one  to  improve  the  inequality 

(2.13)  for  matrices  (i.e.  B  is  finite  dimensional). 


r 


-3  - 


Theorem  2.  Consider  all  polynomials  p (z)  and  q(z)  of  degrees  m  and  n  respectively 
■  ~ 1  trh  that  the  *  unction  f(z)  *  p(z)/q(z)  satisfies  (2.5).  Suppose  that  a  >  1.  Then 
there  exitita  a  positive  constant  p(a,m,n)  such  that 
C.M)  la  |  «  Kpta.m.nlv*"'1 ’  . 

V 

To  prove  this  theorem  we  need  the  following  lemma. 

Lemma  2.  Let  p(z)  be  a  polynomial  of  degree  m.  Then  there  exists  a  constant  K( m) 
such  that 

(2.15)  max  |p(z)|  <  K(m)(  max  I p( re^ ^ > I >  . 

I  z  |  =r  |  6 1  <tt/4 

Proof .  It  is  enough  to  consider  the  case  r  *  1  with  p(z)  of  the  form 

m 

p(z)  -  I  |  (z  -  c.)»  kj  <  Ic.l  •••  *  1;  I  • 

i-i  1  1  2  ™ 

For  m  =  1  it  suffices  to  choose  K ( 1  >  -  5.  Let  m  >  1.  Define 

10 

K'(m)  =  max  !  max  |p(z)|/  max  lp(e  )()  • 

o<lt,  |w*1;  K3  |zl»i  1 6 1 < rr 

i  m 

m-1 

In  case  that  I;  I  >  1  let  q(z)  -  1  f  (z  -  Then 

"  1-1  1 

max  lp(z)l  >  (It  I  +  1'  max  lq(z)l 
lzl-1  m  I z I  —  1 

s  2( I t^l -1 )K(m  -  1)  max  lq(e  )|  4  2K(m  -  1)  max  I p( ) 1  . 

’  ICM‘x/4  1 0 1 F  x  /4 


Put 


K(m)  -  max(KMm),  2K(m  -  1)) 
and  the  lemma  follows.  □ 

Proof  of  Theorem  1 .  Without  restriction  in  generality  we  may  assume  that  p(z)  and 
7(z)  do  not  have  common  zeros.  Also  it  is  enouqh  to  consider  the  case  K  *  1.  The 
inequality  (2.5)  implies  that  we  can  choose  q  and  p  of  the  form 


-4  - 


(2.16) 


p(  z) 


£  — 

Z  A  1  T  < 1  -  zu  ) ,  q(z)  »  1  T  (1  -  Z?,)  . 
i-1  1  i-1  1 


The  inequality  (2.5)  yields  |{  |  <  1,  i  ■  1,...,n.  Put 


m~£  n 

(2.17)  g(z)  =  A  TT  (*  -  u.  )/TT  (z  -  t.  )  , 

i-1  1  i-1  1 


(2.18) 

Also 


lg(z)|  <  |z|m  n+a/(|z|  -  1)a,  |z|  >  1  . 


(2.19) 


g(z) 


00 

I 

v-0 


-(v+n-m) 

a  z 


Izl  >  1  . 


Note  that 


(  2tt  i ) 


-1 


,  (v+n-m-1) 

g(z)z  dz 


■R>1 


Let  D.,...,D  be  p-mutually  disjoint,  open  and  bounded  domains  with  the  boundary 


T . r  respectively.  Assume  that  C.  e  I  1  Dj(  i  -  1,... 

1  p  1  fci  i 


n.  Then  we  obtain 


(2.20) 


l  (2Z1)-1  /  g(z)z(V+n-m-ndz  . 
j-1  r . 


To  obtain  the  estimate  (2.14)  we  are  going  to  choose  the  domains  D1r...,Dp  according  to 
the  configuration  of  and  the  value  of  v.  First  we  group  the  points 

1  ft 

S .,..., S_  following  Horton  (1964).  Let  be  one  of  the  points  with  the  largest 

*1 

modulus,  1 5^  I  -  1  -  S.|  >  I I  <  i  “  1,...,n.  Then  we  form  S1  from  all  those  points 

which  can  be  joined  to  ?  by  a  chain  of  points,  each  link  of  which  has  length  <  {  . 

11  1 
In  the  same  way  S2  is  formed  from  the  remaining  points,  and  so  on  until  all  the 

points  have  been  included  in  some  S„.  For  each  s„  we  denote  by  1-6.  and  i  - 

d  d  cr 

the  modulus  of  the  largest  and  the  smallest  Ic. I,  i.  t  S-.  We  rename  so  that 

lie  i  n 


-5- 


(2.21  )  0  <«<•••<  5  . 

1  s 

Consider  any  particular  s.  and  let  us  denote  its  members  by  X^,  i  =  1,2,...,k,  where 
1  -  j  <  |X. |  v  1  -  5^,  i  =  Let  us  also  denote  the  points  not  in  Sg  by  u  , 

j=1,2,...,n-k.  We  claim 

(2.22)  6.  <  1  -  HI  <  k5  ,  IX.  -  x.l  <  (k  -  1)5,  |X.  -  u.|  >  50  . 

b  i  3  13 

Indeed,  the  first  two  inequalities  follow  immediately  from  the  assumption  that  there  exists 

a  chain  of  at  most  k  points  between  X.  and  X.  such  that  the  distance  of  any  link 

1  3 

<  5,.  The  last  inequality  is  a  consequence  of  u.  not  being  in  s  .  Let 
D  j  S 

m-£ 

(2.23)  h(z)  =  A  |  |  (z  -  u. )  . 

i  =  1  1 

For  X.  f  S  put 
t  D 

(2.24)  n  =  (1  +  25  ) X  / | X  |  . 

D  t  t 

Then 

m-£ 

(2.25)  h(z)  *  l  h . ( z  -  n)3  . 

We  now  estimate  hs.  Let  T  be  a  circle  |z  -  n|  s  Then 

J  D 

Iz  -  ;.|  <  I z  -  *tf  +  Ut  -  c.l  <  \z  -  nl  +  In  -  Xfc|  +  |Xfc  - 
=  5g  +  1  +  26g  -  |Xtl  ♦  |Xt  -  C.l  <  (k  +  3 )60  +  -  ;.l 

where  the  last  inequality  follows  from  (2.22).  In  particular 

!z  -  X.l  <  2(k  +  1)6 
J  “ 

in  view  of  (2.22).  Apply  the  Cauchy  formula  for  hj  and  use  (2.18)  to  get 


-6- 


Ih.l  =  (21T)'1  |/  h(z)dz(2-n)_(  j+1>  I  4  6‘,^+a)4,trt'a  TT  [  (k+3  )6.  +  I  X  -C  || 

3  r  s  i-i  B  t  1 

<  [2  (k  +  i)]k4m+a5"(j+a>+k  TT  [  ( k  +  3)6.  +  IX  -  u .  I) 

B  1,1  B  t  i 

4  [2  (k  +  2)]n4m+Q6-^+a,+kfl  IX  -  m.I  . 

P  .  t  1 


We  now  consider  the  following  three  cases 
(i)  6  >  1/(2n+2nv)  , 

P 

(ii)  6^  4  1/(4nv)  , 

(iii)  neither  (i)  nor  (ii)  holds. 

Here  v  is  a  positive  integer  and  v  >  m  -  n  +  a. 

Case  (i).  Let  C..  be  a  disc  |z  -  S.l  <  5„/2,  for  C.  e  S„.  Then 
-  i  id  id 


(2.27) 


D  “  U  C  =  U  Ds 

i=l  j-1  3 


where  each  0.  contains  a  subset  of  some  S„  and  D,  T  n  =  *  for  j  *  k.  Let  T  be 
3  B  3  k  ■) 

the  boundary  of  Dj.  Then  i(f J  -  the  length  of  r  -  satisfies  the  inequality 


(2.28) 


MT.)  <  2nn(P.)6,  , 
3  3  8 


where  n(D.)  is  the  number  of  points  C,,...,C  in  Dj .  Let  z  e  T..  clearly 
j  1  n  j  j 

z  =  Xt  +  P,  Ip!  -  <53/2#  Sg  =  By  the  definition  of  D  j ,  |z  -  X^|  >  6^/2  for 

1  <  j  *  k.  Also 


|z  -  u.l  -  U  -  u.  +  p|  >  U  -  p.|  -  5/2  >  1  \\  -  VJ .  1  • 

3  tj  t  j  &  2  t  j 


Thus 


(2,29) 


n  „  _v  n-k 

ITT  (z  -  c.)  I  4  2n«  k  TT  lx  -  u, 

i«i  1  j=i  k  3 


,  -1 


Also  for  r;  of  the  form  (2.24)  we  have 


-7  - 


U  -  nl  <  |z  -  X  |  +  IX  -  nl  <  —■  +  1  +  26  -  | \  |  <  (k  +  3)6  . 

t  t  £  p  t  p 

Combine  (2.25 )-( 2 .26 )  with  the  above  equality  to  deduce 
( 2 . 3u  ;  |h(z)l  <  [2(k  +  2)ln+mm  4m+°6^"a  TT  U  -  u  I  . 

8  i-1  1 

Finally  we  deduce 

(2.31)  lg(z)  I  <  (16(n  +  2))n+m+a6'a  . 

P 

Uslnq  the  equality  (...20)  and  the  inequalities  (2.28),  (2.31)  for  v  >  m  -  n  we 

la  |  <  ^  (2,)-1  /  |g(z)||z|  ( V+W':n' U  J  dz  ]  <  n;i6(n  2)Jn+B+a(  min  6.) 

v  j-i  r  i<e<s  8 

<  na(16(n  +  2  „  n+^2  < n+2  )  ( a-1 )  v«-1  # 


as  a  >  1. 
Case  (ii). 
by  (2.27). 
(2.32) 


Thus  we  have  shown  (2.14)  (K  =  1). 

Let  be  an  open  disc  with  center  at  and  radius  1/2v. 

Assume  that  z  e  T  .  So 

z  =  ci/l(;ii  +  P,  IpI  -  ~  ■ 


We  now  estimate 

K( T)  =  max  |h( z) |  , 
zeV 


_1_ 

2v 


i0. 
e  )  , 


191  (5  . 

4 


According  to  (2,18) 


n 

K(D  <  e(4v)a[max  1  I  |z  -  ;  |]  , 

zeT  t=1 


for  v>  >  m  -  n  *  a.  Let  *  M  ♦  —  Clearly  e 

z  e  T  or  z  of  the  form  (2.32)  which  is  in  we  have 


r  . 

j 


We  claim  that 


It. 

l 


< 


Indeed  it  is  easy  to  see  that  for  such 


iz  -  ctl  <  3ln.  -  ctl  . 
z  the  following  inequalities  hold 


get 


a+1 


Form  r> 


for 


-8  - 


So 


lz  -  't'  *  h'  !ni  ‘  V  *  Tv*  ,z  -  v 


i 

<  - 


ln1  -  Cfcl  <  |z-nil+|z-Ctl<4!z-Ctl+lz-;tl=5|z-;tl  , 

|z  -  Ctl  <  \z  -  nil  +  In.  -  Ctl  <  2!  n4  -  it  I  +  In.  -  ;t  I  =  3|n.  -  Cfc  I  . 

Therefore 


KIP)  4  e3n<4v)a  |  I  h  -  ;  I  . 

t-1  1  * 


Let  z  =  CjL/IC±l  +  D  f  I\  ,  |  d|  *  — .  Then  by  Lemma  2  and  the  above  inequalities 


|g(z)|  <  K(nK(?  ~  ^  <  K  ( m  -  t )  ( 1 5  )  n  ( 4v )  °e 
n 


t-i 


|Z  -  ;tl 


and 


|g(z)z''*m_n'1|  <  K(m  -  t )  1 1 5 )  " (4  v)  3e  ( 1  +  J_,  v*ra-n-1  t  K(m  _  £  j  ( 15 ,  n(4  j)  3e~ 
for  v  >  m  -  n  +  a.  As  the  lencrth  of  the  boundary  of  D  does  not  exceed  fro 

(2.20)  we  get 


Case  ( iii  )  . 
(2.33) 


la  I  <  K f ip  -  tlndSl^V^'’ 


In  this  case  we  claim  that  there  exists  1  <  '*  <  s 
1 


6  < 

3+1  „n+2 


vn 


max  c . , 
S  D 


a  n  , 


1 


such  that 


and 

(2.34) 


Y+1 


max  c  .  , 


otherwise  either  (i)  or  (ii)  hold.  (Note  the  inequality  2.21).  Put 

1 

f  2 . 3S )  r  =  (  max  €  .)  +  — 7c —  . 

„  1  n+3 

(K]<y  '  2  vn 


v 


It  is  not  difficult  to  show  that  r  <  Let  5.  r  s  r0r  $  n  v  denote  by  C.  a  disc 

2  V  1  p  1 

with  center  at  and  radius  r.  For  8  >  >  let  be  a  disc  with  center  at 

and  radius  5^/2.  As  before  define  D  by  (2.27).  Now  estimate  ay  from  the  equality 
(2.20)  using  the  arguments  of  the  Cases  (i)  and  (ii)  in  accordance  with  B  >  y  or  B  <  y 
to  deduce  (2.14).  This  concludes  the  proof  of  Theorem  2.  Cl 

Remark  1.  A  special  case  of  Theorem  2,  namely  a  «  1  and  m  =  n  -  1  was  established  in 
Morton  (1964)  . 

Proof  of_  Theorem  1  . 

(A)  =**=>  (P).  Follows  immediately  from  (2.3). 

(K)  ==>  (A).  Let  (I  -  ZA)*1  =  (f  (*))". 

i]  1 

Then  f  A  j <  z  1  =  Pt  ^  (  2) /q1  ^ ( 2l ,  where  the  degrees  of  and  are  n  -  1  and  n 

respectively.  Now  (1.3)  follows  from  Theorem  2. 
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